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1. Introduction 



In this paper, we are investigating the existence of isolated excited states in certain one-dimensional, 
quantum spin models of magnetic systems. It turns out that if the spins are of magnitude 1 or more and 
their interactions have a suitable anisotropy, such as in the ferromagnetic XXZ Heisenberg model, isolated 
excited states are possible. For the spin 1/2 chain the ground states are separated by a gap to the rest of 
the spectrum, and there are no isolated eigenvalues below the continuum. 

Our main result is a mathematical demonstration that such states indeed exist for sufficiently large 
anisotropy. Concretely, we study the one-dimensional spin J ferromagnetic XXZ model with the following 
boundary terms. The Hamiltonian is 



where S^,S^ and are the spin J matrices acting on the site a. Apart from the magnitude of the spins, 
J, the main parameter of the model is the anisotropy A > 1, and we will refer to the limit A — > oo as 
the Ising limit. In the case of J = 1/2 these boundary conditions were first introduced in [2]. They lead 
to ground states with a domain wall between down spins on the left portion of the chain and up spins on 
the right. The domain wall is exponentially localized [T]. The third component of the magnetization, M, 
is conserved, and there is exactly one ground state for each value of M. Different values of M correspond 
to different positions of the domain walls, which in one dimension are sometimes referred to as kinks. In 
[3] and [4] the ground states for this type boundary conditions were further analyzed and generalized to 
higher spin. A careful analysis of the Ising limit, see Section [4l reveals that for J > 1 one or more low- lying 
excitations, each with a finite degeneracy, closely resemble the domain wall, i.e. kinked, ground states, and 
therefore one should expect them to be resolute under perturbations. In this paper, we first show that these 
states exists and correspond to isolated eigenvalues of the finite volume XXZ chain with sufficiently strong 
anisotropy. We illustrate this feature in Figure [1] Moreover, as consequence of the strong localization near 
the position of the ground state kink, these eigenvalues only weakly depend on the distance of the domain 
wall to the edges of the chain, and for this reason, we are next able to demonstrate that they persits even 
after the thermodynamic limit. The main difficulty we must overcome corresponds to the fact that, in the 
thermodynamic limit, the perturbation of the entire chain is an unbounded operator, and therefore, the 
standard, finite-order perturbation theory is inadequate for a rigorous argument. 

The XXZ kink Hamiltonian commutes with the operator Sf ot = ^2a=—L ^a- We define Hm to be the 
eigenspace of Sf ot with eigenvalue M e {— J(2L + 1), . . . , J(2L + 1)}. These subspaces are called "sectors", 
and they are invariant subspaces for i/^(A _1 ). 

It was shown in [3j [5J [6] that for each sector there is a unique ground state of H^(A~ 1 ) with eigenvalue 
0. Moreover, this ground state, ipM, is given by the following expression: 



where the sum is over all configurations for which ^ Q m a = M and the relationship between A > 1 and 
q G (0,1) is given by A = (q + q~ 1 )/2. A straightforward calculation shows a sharp transition in the 
magnetization from fully polarized down at the left to fully polarized up at the right. For this reason they 
are called kink ground states. In [5] Koma and Nachtergaele proved that the kink ground states (as well as 
their spin-flipped or reflected versions the antikinks) comprise the entire set of ground states for the infinite- 
volume model, aside from the 2 other ground states: the translation invariant maximally magnetized and 
minimally magnetized all +J and all -J groundstates. Since the infinite-volume Hamiltonian incorporates all 
possible limits for all possible boundary conditions, this is a strong a posteriori justification for choosing the 
kink boundary conditions. It is worth noting that it has been proved for the antiferromagnetic model that 
no such ground states exist [7J [5] . 

In [9], Koma, Nachtergaele, and Starr showed that there is a spectral gap above each of the ground states 
in this model for all values of J. Based on numerical evidence, they also made a conjecture that for J > | 
the first excited state of the XXZ model is an isolated eigenvalue, and that the magnitude of the spectral gap 
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Figure 1 . The ground state and first excited state of the XXZ chain of length 7 with J = | , 
A = 2.5 in the sector M = -3/2. 



is asymptotically given by Jj(A), where 7(A) is an eigenvalue of a particular one-particle problem. Caputo 
and Martinelli [10] showed that the gap is indeed of order J. 

Our main result is a proof that for all sufficiently large A the first few excitations of the XXZ model are 
isolated eigenvalues. This is true for all J > | and for spin 1 with M even, which is illustrated in Figure [1] 
for the spin 3/2 and spin 2 chains. See Section [3] for the precise statements. It turns out that in the Ising 
limit the eigenvalues less than 2 J are all of multiplicity at most 2 in each sector. Moreover, the first excited 
states are simple except in the case when J > 1 is an integer and M = mod 2J. In this case, they are 
doubly degenerate. This is discussed in Section^ In Section[51 we write the XXZ Hamiltonian as an explicit 
perturbation of the Ising limit. Theorem 1 5 . 31 verifies that the perturbation is relatively bounded with respect 
to the Ising limit, and we finish this section by demonstrating that our estimates suffice to guarantee analytic 
continuation of the limiting eigenvalues. It is clear that the same method of proof can be applied to other 
Hamiltonians. 

While the question of low-lying excitations is generally interesting, it is of particular importance in the 
context of quantum computation. For quantum computers to become a reality we need to find or build 
physical systems that faithfully implement the quantum gates used in the algorithms of quantum computation 
jllj . The basic requirement is that the experimenter has access to two states of a quantum system that can be 
effectively decoupled from environmental noise for a sufficiently long time, and that transitions between these 
two states can be controlled to simulate a number of elementary quantum gates (unitary transformations). 
Systems that have been investigated intensively are single photon systems, cavity QED, nuclear spins (using 
NMR in suitable molecules), atomic levels in ion traps, and Joscphson rings [T^l H31 HI] ■ We believe that 
if one could build one-dimensional spin J systems with J > 3/2, which interact through an anisotropic 
interaction such as in the XXZ model, this would be a good starting point to encode qubits and unitary 
gates. The natural candidates for control parameters in such systems would be the components of a localized 
magnetic field. From the experimental point of view this is certainly a challenging problem. This work is a 
first step toward developing a mathematical model useful in the study of optimal control for these systems 
such as has already been carried out for nuclear magnetic resonance (NMR) [15] [16] and superconducting 
Joshepson qubits [17] , 



2. Set-up 

We study the spin J ferromagnetic XXZ model on the one-dimensional lattice Z. The local Hilbert space 
for a single site a is TL a = C 2J+1 with J e = {0, |, 1, |, 2, . . . }. We consider the Hilbert space for a finite 
chain on the sites [-L, L] = {—L, —L + 1, . . . , +L}. This is H[-l,l] = ®a=-L The Hamiltonian of the 
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Figure 2. The low- lying spectrum of the XXZ model for J — | and J = 2 for various 
sectors. Note the two-fold degeneracy for spin 2 in the sector M = corresponding to 
E + (0, 1) = E_(0, 1). The shaded region contains those eigenvalues which converge to con- 
tinuous spectrum in the thermodynamic limit. We note that, for A = 1, the gap vanishes 
in the thermodynamics limit. 



spin- J XXZ model is 

L-l 

H L (A-') = £ Va+itA" 1 ), 

where S^S^ and S 1 ^ are the spin-J matrices acting on the site a, tensored with the identity operator acting 
on the other sites. The main parameter of the model is the anisotropy A > 1 and we get the Ising limit as 
A — > oo. It is mathematically more convenient to work with the parameter A -1 , which we then assume is 
in the interval [0, 1]. As we said, A -1 = is the Ising limit, and A -1 = 1 is the isotropic XXX Heisenberg 
model. It was shown [2j El HJ E] that additional ground states emerge when we add particular boundary 
terms. Examples of this are the kink and antikink Hamiltonians 



(2.2) HliA- 1 ) = H^A- 1 ) + Jy/1 - A-2 (S 3 _ L - S 3 L j 



(2.3) H? (A- 1 ) = H L {A~ X ) - Jy/l - A-2 (S 3 _ L - S 3 L ) 

It is easy to see that the kink and anti-kink Hamiltonians are unitarily equivalent. We will be mainly 
interested in the kink Hamiltonian with L > 2. Note that, by a telescoping sum, we can absorb the 
boundary fields into the local interactions: 

(2-4) H^A- 1 ) = J2 hl a+1 (A-i), 

a = -L 



The Ising kink Hamiltonian is the result of taking A — > oo, equivalently setting A -1 = 0, namely it is 



L-l 

a=-L 



(2-5) 

= (J + S 3 a )(J - S 3 a+1 ) . 
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Each of the Hamiltonians introduced above commutes with the total magnetization 

L 

■Stat = J]] Sa ■ 
a=-L 

As indicated in the introduction, for each M € {— J(2L + 1), . . . , J{2L + 1)}, the corresponding sector is 
defined to be the eigenspace of Sf ot with eigenvalue M; clearly, these are invariant subspaces for all the 
Hamiltonians introduced above. 

The Ising basis is a natural orthonormal basis for TL[-l,l]- At each site we have an orthonormal basis of 
the Hilbert space TC a given by the eigenvectors of and labeled according to their eigenvalues. We will 
denote this by S^\m) a = m\m) a for m G [— J, J], and a G [—L,L\. Here, and throughout the remainder of 
the paper, we will use the notation [— J, J] for the set {—J, — J + 1, . . . , J} as we have done with [-L, L]. 
Finally, there is an orthonormal basis of the entire Hilbert space consisting of simple tensor product vectors: 

<S>a=-L \ m *)<±- 

Also recall that the raising and lowering operators are defined such that 

+ j \J J{J + 1) - m{m + 1) |m + 1} Q if - J < m < J - 1, 

o a \m) a — < 

10 it m = J; 



f ^JJ{J + 1) - m(m - 1) |m - l) a if-J + 



1 < m < J, 
if m = - J. 



A short calculation shows that S 1 ^ and S 1 ^ are given by S\ = (5+ + S~)/2 and = (S*+ — S~)/2i, and 

^a^a+l + S a S a+1 = — (S^S a + 1 + S a S a+1 ) . 

3. Main theorem 

Many of the results in this paper concern the kink Hamiltonian given by (|2.4[) . We will study it as a 
perturbation of the Ising Hamiltonian (|2.5p in the regime < A -1 -c 1. We denote an Ising configuration 
as m = (m a )^ = _ L , where m a G [—J, J] for each a, and the corresponding basis vector as 



a=-L 

Observe that the Ising kink Hamiltonian is diagonal with respect to this basis, 

#l(0) \m) = E k (m) \m) , where 

L-l 

(3.1) E k (m) = ^2 e k (m a ,m a+1 ) and 

a— — L 

e k (m Q ,TO Q+ i) = (J + m a )(J - m a+1 ) . 

Since each of the e k (m Q , m a +i) are non-negative, it is easy to see that the ground states of H^(0) are all of 
the form 

for — L < a < x — 1 

(3.2) ^o(x,m; L) — \rh) , where m a = ^ m for a = x 

for x + 1 < a < L, 

with some x G \—L,L] and m G [—J, J]. Note that the total magnetization corresponding to ^o(x,m; L) is 
M = —2Jx + m. As we will verify in Proposition [47T1 it is easy to check that these ground states are unique 
per sector. We do point out that there is a slight ambiguity in the above labeling scheme, however, since 
^oi x ~ 1) ~ J'i L) and ^('-'i J\ L) coincide. Let us consider the following elementary result. 

Theorem 3.1. For J G ±N, L > 2 and M G {~J{2L + 1),...,J(2L + 1)}, the eigenspace of if£(0) 
corresponding to energy E = 2J has dimension at least equal to 2L + 1. 
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Figure 3. The two kinds of low lying Ising excitations. The left picture belongs to the set 
K+ (m) and the right to the set K- (m) . 



Because of this theorem, we see that, in the L — » oo limit, the energy E = 2J is essential spectrum; an 
eigenvalue with infinite multiplicity. As our interest is in perturbation theory, it is natural for us to restrict 
our attention to energies strictly between and 2 J. We will call the corresponding eigenvectors "low energy 
excitations" . 

We will now describe all the low energy excitations for the Ising kink Hamiltonian. In order to do so, it 
is convenient to introduce a family of eigenvectors which contains all possible low energy excitations. 

Let x be any site away from the boundary, i.e. take x £ [— L + 1,L — 1], and choose to £ [—J, J]. These 
choices specify a sector M = —2Jx + to and a groundstate 1 ^>o(x,m; L). If to < J, we define the following 
vectors: For 1 < n < J — to set 



(3.3) ^>^(x,m;L) = \rh) , where m a - 
If — J < to, we define the following vectors: For 1 < n < J 

(3.4) ty~(x, to; L) = \rh) , where m a = 



We call these two sequences of Ising basis vectors the "localized kink excitations" . Clearly, these vectors 
{^^(x, to; L)} have the same total magnetization M and moreover, 

(3.5) Hl(0)^(x,m;L) = E±(m, n) ^„(x, m; L), with E±(m, n) = n 2 + (J±m)n, 

where E±(m,n) has been calculated using ()3. 1|) : note that in each case there is only one non-zero term. In 
Figure[3]we have shown two ground states, in the left and right graphs, as well as two low excitations. These 
are the schematic diagrams for ^o{ x : m )- m solid in both pictures, and $>t {x, m) and 'J ^7 (x, m), in a dashed 
line, in the left and right, respectively 
We now define two sets of labels: 

(3.6) K±(m) = {n : n£N, 1 < n < J^m, E±(m,n) < 2J} . 

Depending on to and J, neither, one, or both of these sets may be empty. We have the following theorem. 

Theorem 3.2. (1) The low energy excitations of H\(0) form a subset of the localized kink excitations 
introduced in Ii3.3\) and j3.4\ l above. 

(2) For M of the form M = —2Jx + to with some x £ [-L + 1, L — 1] and to £ [—J, J], the set of low energy 
excitations equals 

{*+(x,m;L) : n £ K + (m)} U {^~(x, to; L) : n £ K_(m)} . 
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This is a nonempty set except in the following two cases: J = ■k, or J = 1, and to = 0. 

(3) The low energy excitations of H^(0) are at most two-fold degenerate. The first excitation is simple, 

except for the case that J is an integer J > I, and M = mod 2 J . 

Remark 3.3. The two-fold degeneracy of the first excited state, guaranteed by part (3) of the theorem, 
occurs due to the spin flip and reflection symmetry. All other degeneracies with energy < 2J occur as 
follows. Suppose 2 J = ab for integers 2 < a < b. In this case, let m = J — 2 + a — b. Then "^(x, to; L) is 
degenerate with ^I'~_ 1 (.t, to; L), both with energy 2 J — 1 + a — b. Similarly, ^i(x, —m; L) is degenerate with 
^'^_ 1 (x, — m; L) for the same energy. 

Next we consider the perturbed Hamiltonian iJ|(A _1 ). As is discussed in Section T5.21 each low- lying 
eigenvalue E±(m,n), associated to some n G K±(m), is isolated from the rest of the spectrum by an 
isolation distance d±{m, n) > 1, independent of L, which is defined by 

d±(m,n) = inf \E — E±(m, n)\ . 

EE<r(H*(0))\{E±(m,n)} 

Theorem 3.4. Let L > 2, and fix J > |. For any n G K±(m), consider the interval I = \E±(m,n) — 
d±(m,n)/2, E±(m,n)+d±(m,n)/2] about the low lying energy E±(m,n). The spectral projection of H k (A^ 1 ) 
onto I is analytic for large enough values of A. In particular, the dimension of the spectral projection onto 

1 is constant for this range of A. 

Remark 3.5. Our estimates yield a lower bound on A as is provided by \5. 39$ . A slightly worse bound 
demonstrates that taking A > 18 J 5 / 2 suffices, but we do not expect either estimate to be sharp. 

The above theorem confirms the structure of the spectrum shown in Figure [1] Moreover, since our 
numerical calculations indicate that some of the eigenvalues enter the continous spectrum, we do not expect 
this type of perturbation theory to work for the entire range of A -1 G [0, 1]. 

4. Proof of Theorems 13.11 and 13.21 (Excitations of the Ising Model) 

In this section, we will focus on the Ising kink Hamiltonian i?£(0), as introduced in (|2.5p . for a fixed 
J G iN. The ground states can be characterized as follows: 

Proposition 4.1. The ground states of the Ising kink Hamiltonian are all of the form 

{—J for a = —L, . . . , x — 1, 
m for a = x, 
+J for a = x + 1, . . . , L; 

for x G [-L, . . . , L] and to G [—J, J]. Moreover, there is exactly one ground state in each sector; the total 
magnetization eigenvalue for 'I'ot^j to; L) is M = 2Jx + m. 

This proposition was already used, implicitly, in setting up Theorem 13.21 Here we prove it. 
Proof: Given any Ising configuration to = (m a )a=-L' equation (|3.1[) demonstrates that the energy 
associated to |m), E k (m), is the sum of 2L non-negative terms. Therefore, the only way to have E k (m) = 
is if all of the summands are 0. This is the case only if either m a = — J or m a +i = +J for all a. Clearly, 
this is satisfied for the Ising configurations with m a = — J for all a < x, m a = +J for all a > x, and m x 
equal to any number in [—J, J]. It is equally easy to see that these are all of the ground state configurations: 
if to is a groundstatc configuration and for some x we have m x =/= — J, then m x+ \ = J, which, by induction, 
means that m a = +J for all a > x. Similar reasoning yields that m a = — J for all a < x. 

To show that the ground states are unique in each sector, consider the equation M — to = —2Jx, subject 
to the constraints: M G {- J(2L + 1), . . . , J(2L + 1)}, to G [- J, J], and x G [-L, L]. If M - J + mod 

2 J, then there is a unique pair (x,m) which satisfies this equation. If M — J = mod 2 J, then there are 
two possible solutions (x, J) and (x — 1, — J) for some x. But then, it is trivial to see that ^0(^7 J] L) and 
^o(x — 1, —J; L) coincide. | 

Thus in any sector M, there is a unique ground state eigenvector ^o(x,m\L) for some choice of (x,to) 
with M = —2Jx + m. We next observe that there arc many eigenvectors with eigenvalue 2 J; recall this was 
the statement of Theorem 13.11 
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Proof of Theorem I3.lt First consider the case that M — J is not divisible by 2 J. Then the unique 
groundstate eigenvector is *&o(x,m; L) for some x G [—L,L] and — J < m < J. If x > — L, then for each 
y G [— L, x — 2] consider the Ising configuration to with components 

— J if a < x — 1 and a ^ y; 

— J + 1 if a = y; 
to — 1 if a = x; 
J ifa>x+l. 

The total magnetization for the vector \rh) is still M. Using the formula (|3.1|) again, it is easy to check that 
e k (rn y , TOy+i) = 2 J, e k (m a , m a +i) = for all a ^ y, and therefore, E k (m) = 2 J. This constitutes x — 1 + L 
possible values of y; producing at least this many distinct eigenvectors with energy 2 J. Similarly, if x < L 
then there arc L — x — 1 Ising configurations of the form to with components 




if a > x + 1 and a ^ y; 

corresponding to some y G [x + 2,L]. In total, this gives 2L — 2 orthonormal eigenvectors corresponding 
to eigenvalue 2J, yielding a lower bound on the dimension of the eigenspace. If x = —L or L, then the 
dimension is increased by at least 1. In the special case where M — J is divisible by 2 J, the eigenvector can 
be written in two ways as ^>o{x, J; L) or ^o(x — 1, — J; L). When constructing excitations for y to the left of 
the kink, use the first formula above relative to *&o(x, J; L). When constructing excitations for y to the right 
of the kink, use the second formula above relative to ^o(x — 1, — J; L). Once again, this results in 2L — 1 
orthonormal eigenvectors corresponding to eigenvalue 2 J. | 

We now claim that any Ising configuration which is neither a ground state nor a localized kink excitation 
corresponds to an energy that is at least 2 J. This is the content of the following lemma. 

Lemma 4.2. Consider an Ising configuration rh = (m a )^ = _ L . 

(1) If there is any x G [— L, L — 1] such that m x > to^+i, then E^{m) > 2 J. 

(2) If there is any x G [— L + 1, . . . , L — 1] such that m x -\ > — J and to^+i < J, then E k (m) > 2 J. 
Proof: (1) It is clear from (|3.ip that we need only prove 

(4.2) J 2 - m x m x+ i + J(m x - m x+ i) = e k (m x ,m x+ i) > 2.1. 

Since m x > m x+ i, we have that J(m x — m x+ i) > J. We need only verify that J 2 — to^to^+i > J to establish 
the claim. The product of two integers m x , m x+ i G [—J, J] is at most equal to J 2 . But this is only attained 
by m x = m x +i = ±J. Since m x > m x+ i, we can neither have m x +\ = J nor m x = —J. The next largest 
possible value of to^to^+i is J( J — 1), and we have verified the claim. 
(2) Again, because all terms are nonnegative, it suffices to show that 

e k (m x ^i, m x ) + e k (m x , m x +i) > 2 J . 

Using the formula for e k (-, •) and simplifying gives 

e k (m x _i,m x ) + e k (m x , m I+1 ) = 2J 2 + J(m x ^i - m x+ i) ~ m x (m x _ 1 + m x+1 ) . 

Since \m x \ < J, we have 

e k (m x -i, m x ) + e k (m x , m x +x) > 2 J 2 + J(to. t _i - to^+i) — J|to x _i + m x+ i\ 



Therefore, 



= J(2J + m x -i - m x+ i - |m x _i + m x +i|) 
e k (m x -i,m x ) + e k (m x ,m x+ i) > 



2 J( J — m^x) if to x _i + m x+ i > 0; 
2J(J + to t _i) if to x _i + m x+ i < 0. 
Since J — m x+ \ > 1 and J + m x ~\ > 1, in either case we have proven the claim. 
Now we can finish the proof of the first main result, Theorem 13.21 
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Proof of Theorem 13. 2t Part (1) of the theorem is a direct consequence of Lemma T4.2I because the only 
Ising configurations that do not satisfy condition (1) or (2) of that lemma are the ground state configurations 
and the localized kink excitations. 

To prove part (2), let M = —2Jx + m for some x G [— L + 1, L — 1] and m G [—J, J]. We first consider 
the case that J > 3/2. 

If J is sufficiently large and \m\ < J — 2, then n = 1 G K±(m) as 

£±(m,l) = l 2 + (J±m)l < 1 + J+H < 2J. 

So, in this case, both K+(m) and K-(m) are nonempty. 
Similarly, if |m| = J — 1, then cither 

E+(-J + 1,1) = l 2 + (J+(-J + l))l = 2, 



or 

E-(J - 1, 1) = l 2 + (J - (J - 1))1 = 2. 

Hence, 1 G K + (m) U 7\_(m) if \m\ = J - 1 and J > 3/2. 
Lastly, if |m| = J, then either 

£+(-J,l) = l 2 + (J+(-J))l = 1, 

or 

£_(J,1) = l 2 + (J — J)l = 1. 

Hence, 1 G K+(m) U K_(m) if \m\ = J - 1 and J > 1. 

We have proved (2) in the case that J > 3/2. Actually, the last observation above also verifies (2) in the 
case that J = 1 and m = ±1. 

Finally, if J = 1 and m = 0, then E±(0, 1) = l 2 + (1 ± 0) = 2 and if J = 1/2, then E±(+l/2, 1) = 
l 2 + (1/2 — 1/2) = 1. In both these cases, the set of kink excitations is empty. 

We now prove the first part of (3). First, observe that for any m G [— J, J], J± m > and therefore with 
m fixed, both E±(m,n) are increasing functions of n for n > 0. Therefore the only degeneracies that can 
occur for a particular energy E is if E+(m,ni) = E and £L(m,n2) = i? for some integers n\ and ri2. This 
is obviously at most a two-fold degeneracy. 

In order to prove the second part of (3) , we will simply prove Remark 13.31 Without loss of generality, 
suppose m > 0. Then E+(m,2) = 4 + 2 J + 2m > 2 J. So the only possibility for E+(m,n) to be less than 2 J 
is if n = 1, which gives the energy J + m + 1 and an auxiliary condition, namely J — m > 2. A degeneracy 
happens only if there is a p > 1 such that E-(m,p) = E + (m, 1). This means 

p 2 + (J - m)p = J + m+ 1 -i==^> (p+l)(p + J-m-l) = 2J. 

Setting a = p + 1 and b = p + J — m — 1 (which satisfies b > a because J — m > 2) we have exactly the result 
claimed. Note that the second part of (3) refers to the first excitation. For m > 0, the lowest excitation is 
E-(m, 1). This means p = 1, which implies m = and therefore 2 J is even with J > 2 (because J — m > 2). 



5. Proof of the main theorem 

The goal of this section is to prove Theorem l3.4l We will do so by analyzing the kink Hamiltonian if£(A -1 ) 
as a perturbation of the Ising limit H^(0). Within the first subsection below, specifically in Theorem 15. 3[ 
we prove that the operators which arise in our expansion of iJ^(A _1 ) are relatively bounded with respect 
to ij£(0). In the next subsection, we discuss how the explicit bounds on A, those claimed in Theorem 13. 4[ 
follow from relative boundedness and basic perturbation theory. 
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5.1. Relative Boundedness. In this subsection, we will analyze the kink Hamiltonian introduced in (|2.4p . 
Recall that this Hamiltonian is written as 

(5-1) ^(A- 1 ) = J2 h^ a+1 (A-'), 



a=-L 



^.a+iCA" 1 ) = J 2 - SlSl +1 A'VX+i + S 2 a S 2 a+1 ) + Jy/1 A-2 (S 3 a - S 3 a+1 ). 
By adding and subtracting terms of the form J(S 3 t — to the local Hamiltonians, we find that 

(5.2) HliA- 1 ) = Hl(0) + A^H^ + (l - \fl - A~ 2 ^ 

where 

L-l 

(5-3) H™ = - Co+l with C+i = 2(^+1+^+1) 

a— — L 

and 

(5.4) ff[ 2) = J(5£-5i L ). 

In Theorem 15.31 below, we will show that both and Hjp are relatively bounded perturbations of H^{0). 
To prove such estimates, we will use the following lemma several times. 

Lemma 5.1. Let A and B be self-adjoint n x n matrices. If A > and Ker(A) C Ker(B), then there exists 
a constant c > for which, 

(5.5) - cA < B < cA. 

One may take c = where Ai denotes the smallest positive eigenvalue of A. 

Proof: Any vector ip £ C" can be written as ip = ipo + i>i where ipo £ Kei(A) and ipi £ Ker(A) ± . Clearly 
then, Ai||-0i|| 2 < (ip,Atp) and therefore 

(5.6) = |(^,%)| < ||B||||^i|| 2 < 

as claimed. ■ 

For our proof of Theorem [531 wc nn d it useful to introduce the Ising model without boundary conditions 
as an auxiliary Hamiltonian, i.e., 

h-\ 

h l{0) = Yl where h a , a+1 {0) = J 2 - S 3 a S 3 a+1 . 

a— — L 

It is easy to prove the next lemma. 

Lemma 5.2. The Ising model without boundary terms is relatively bounded with respect to the Ising kink 
Hamiltonian. In particular, for any vector ip, 

\\H L m\\ < mtmw +2j 2 iivii- 

Proof: Consider the terms of the Ising kink Hamiltonian: 

(5.7) < Q+ i(0) = (J + Sl)(J - S 3 a+1 ) = h a . a+1 (0) + J(S 3 a - S 3 a+1 ). 
Summing on a then, we find that 

(5.8) H L (0) = H$(0) + J(Sl-St L ), 
and therefore, the bound 

(5-9) \\Zl(0M < \\H k L m\\ +2J 2 |H|, 

is clear for any vector ip. | 
Wc now state the relative boundedness result. 
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Theorem 5.3. The linear term in the perturbation expansion of H^(A 1 ), see I5.2j) . satisfies 
(5.10) 



< VJ 2 + 2J3 ||^(0)^|| + 2JW.J 2 + 2J 3 Hl 
for any vector ip. Moreover, we also have that 



(5.11) 



4 a V 



< 2J 2 ||VII 



Proof: Using Lemma I5T21 it is clear we need only prove that 
(5.12) 



< VJ 2 + 2J3||i? L (0)V||. 



to establish (|5 . lOf) . To this end, Lemma 1 5. II provides an immediate bound on the individual terms of 

these Hamiltonians. In fact, observe that for any fixed a, both /i a , a +i(0) and h^ a+1 arc self-adjoint with 
h a ,a+i(0) > and 

(5.13) 



Ker (/i a , a +i(0)) = {\rh) : m a = m a+1 = ±J } C Ker \h K a ' a+1 
It is also easy to see that, for every a, the first positive eigenvalue of /i Q , a +i(0) is Ai = J, and we have that 



(5.14) 



\h 



(i) 



2 ll^^a+l + ^+l|| — J"' 



An application of Lemma 15.11 yields the operator inequality 

(5.15) -Jh a . a+1 (0) < h ( ^ a+1 < Jh a , a+1 (0), 

valid for any a. 

The norm bound we seek to prove will follow from considering products of these local Hamiltonians. For 
any vector -0, one has that 

L-l 

r(!)„/,l|2 _ /„,, fzrWV^ _ /„/, J 1 ) 1,(1) 



(5.16) 

and 
(5.17) 



\m 



a, a' — — L 



L-l 



||ff L (0)V|| 2 = (V,(^l(0)) 2 ^) = 2 (^fc«,a+i(0)W + i(0)V>. 



a,a' — ~L 



The arguments we provided above apply equally well to the diagonal terms of (|5 . and (|5 . 1 7[) in the 
sense that 

(5.18) 



J 2 (h a , a+ i(0)) 2 < (h$ a+1 ) < J 2 (h a , a m\ 



is also valid for any a. We find a similar bound by considering the terms on the right hand side of (|5.16[) 
and (|5.17p for which \a — a'\ > 1. In this case, each of the operators h atQ:+ i(0) and h^ a+1 commute 
with both of the operators h a > t0l > +i(0) and h^) a , +1 . Moreover, we conclude from (|5.15[) that the operators 

= J/i Q . Q +i(0) ± h^ a+1 are non-negative for every a. Since all the relevant quantities commute, it is 
clear that 



(5.19) 



< - 0?+9 Q , + g a gt>) = J 2 K, a +i{Q)h a >, a > +1 {Q) - h ( ^ a+1 h^) a , +1 . 



Our observations above imply the following bound 

L-2 



(5.20) 



- j 2 \\h l w\\ 2 < (*> + h^ a+2 h^ a+1 )i 



a=-L 



In fact, the terms on the right hand side of (|5.20p for which either a' — a or \a — a'\ > 1 arc non-positive 
by (|5.18p . respectively, (|5. 19[) . In the case that \a — a'\ = 1, the operators /i Q , Q +i(0)/i Q '. Q '+i(0) are non- 
negative (since they commute) and hence we may drop these terms; those terms that remain we group as 
the self-adjoint operators appearing on the right hand side of (|5.20[) above. 
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Our estimate is completed by applying Lcmma l5.1l one more time. Note that for any a £ {—-£•, ■ 
the operator A a = /i aiQ +i(0) + /ia+i,a+2(0) is self-adjoint, non-negative, and 

(5.21) Ker(A a ) = {\m) : m a = m a+ i = m Q+ 2 = ±J} C Ker(5 Q ), 

where the self-adjoint operator B a appearing above is given by 



,L-2} 



(5-22) ±j a — ' t Q , Q +i' t Q+i, Q! +2 

For each a, the first positive eigenvalue of A a is \\{pt) 



ct+l,a+2 a,a+l " 

J, and it is also easy to see that ||B a || < 2 J 4 . 



Thus, term by term Lemma |5. II implies that 

(5.23) (ip,B a i>) < 2J 3 (V,A Q V), 

from which wc conclude that 



(5.24) 



< (J 2 + 2J 3 ) \\H L {Q)^\Y 



as claimed in f|5. 12|) . We have proved (|5.10|) . 

Equation (|5 . follows directly from the easy observation that 



H 



(2) 



is equal to 2 J 2 . 



I 



5.2. Perturbation theory. In Section 01 we verified that, in any given sector, the spectrum of the Ising 
kink Hamiltonian, ifk(O), when restricted to the interval [0, 2 J) consists of only isolated eigenvalues whose 
multiplicity is at most two. In fact, for the sector M = ~2Jx + m these eigenvalues are determined by 

(5.25) E±(m,n) = n 2 + (J ± m)n 

for those values of n £ N with E±(m, n) < 2 J. It is clear from (|5.25p that each of these eigenvalues have an 
isolation distance d±(m, n) > from the rest of the spectrum and that this distance is independent of the 
length scale L. 

For our proof of the relative boundedness result in Theorem 15.31 we expanded the Hamiltonian as 

r(2) 



(5.26) 



HKA- 1 ) = H k L (Q) + A-'Hi 



1 - \J\ - A- 2 ) H 



Using the first resolvent formula, it is easy to see that 



(5.27) 



i 



A if 



(i) 



1 - y/l- A- 



H 



(2: 



R(0 



where we have denoted the resolvent by i?(£) 
small enough so that 



(if£ (0) — ^) 1 , and it is assumed that A 1 has been chosen 



(5.28) 



(A^ffW + (l - v/T^A^) tff) R(t) 



< 1. 



It is clear from sections II. 1.3-4 of [TB] and chapter I of [115] that the spectral projections corresponding to 
ff^(A _1 ) can be written as a power series in A -1 , the coefficients of which being integrals of the resolvent 
over a fixed contour T. Proving an estimate of the form (|5.28p for A large enough, uniformly with respect 
to £ £ r, is sufficient to guarantee analyticity of the spectral projections. We verify such a uniform estimate 
below. 

Let E±(m,n) be an eigenvalue of if£(0) with isolation distance d±(m,n) as specified above. Denote by 
r the circle in the complex plane centered at E±{m,n) with radius d±(m,n)/2. We claim that if 

(5.29) A > 18J 5/2 , 

then (|5.28p is satisfied uniformly for £ £ T. 

Wc proved in Theorem 15.31 that for any vector tp, 



(5.30) 



< VJ 2 + 2J*\\H$ (0)^|| + 2J 2 VJ 2 + 2J 3 |HI 



Applying this bound to vectors ip of the form ip = R(£)<f) yields a norm estimate on if T f£(£), i.e. 



(5.31) 

Moreover, since 
(5.32) 



< 



VJ 2 + 2J3 \\Hl(0)R(£)\\ + 2J 2 VJ 2 + 2J*\\R(0\\) II 



\\H k L (Q)R(0\\ < 1 + ielTOII 
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we have proved that 



(5.33) 



A-'H^Rit) < A-V^ 2 + 2J 3 [1 + (|e| + 2J 2 ) 11^(011 ] • 
Similar arguments, again using Theorem 15. 3[ imply that 
(5.34) 



(5.35) 

and 
(5.36) 



1 - V1-A-2) H^R(0 < (1 - Vl-A-2) 2J 2 ||ij(0|| 
describee 



For £ G r, the circular contour described above, we have that 

1 2 



dist(£,cr(#£(0))) d±(m,n) : 



|CI < E±(m,n) + d±(m,n)/2. 
We derive a bound of the form (|5.28p . uniform for £ 6 T, by ensuring A large enough so that 
(5.37) CiA- 1 + C 2 (l - y/l- A~ 2 ) < 1, 

where 



(5.38) 



C a = v 7 J 2 + 2 J 3 



1 



(2E±(m,ri) 



1 



4J 2 



\ d±(m,n) d±(m,n) 
Explicitly, one finds that the inequality (|5.37|) is satisfied for all 

cl + cl 



and C2 



4J 2 



d±(m, n) 



(5.39) 



A > 



C 2 + 2C 2 - 1 + d - CiC 2 
Equation (|5.29p is a simple sufficient condition for A to satisfy this inequality. This is easy to verify if one 
first replaces 1 - s/l - A~ 2 by A" 2 in ([537]) . 
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